Two types of models are used to simulate contaminant transport in ground water: models that simulate advection, dispersion, and chemical reactions, and models that simulate only advection and chemical reactions. The theoretical difficulties described above are pertinent to models that incorporate dispersion, and therefore, current dispersion models are not suitable for routine use in predicting ground-water-quality changes. Preliminary research applications of models to dispersion problems related to coal mining are presented in Amend and others (1976) and Homsy (1979); however, for many applications, it is unnecessary to include dispersion in the simulation. Models that include only advection and chemical reactions are described below, and problems relating to coal mining are used as examples.
In the East, the major effect on ground-water quality is acid mine drainage. Investigators at Ohio State University (Morth and others 1972, Ricca and Chow 1974) developed a mathematical model to predict daily mine-water discharge and acid loads and tested the model using data from a deep mine in southern Ohio. However, the model requires detailed input data including precipitation and pan-evaporation records, initial soil moisture conditions, initial ground-water storage conditions, oxidation rate parameters, initial acid storage, and flow and acid load coefficients.
In the West, applications of mathematical models to water-quality problems generally have dealt with predicting water-quality changes as a result of water movement through spoils. McWhorter and others (1979) presented two such mathematical models, each of which performs at a different level of complexity but neither of which includes dispersion. The more complex of the two is a simulation model that allows consideration of flow and contaminant transport through both the unsaturated and saturated zones. All chemical reactions are assumed to take place in the unsaturated zone. The following processes are considered in the model: solubility and precipitation of gypsum, chemistry of undissociated calcium and magnesium sulfate, exchange of calcium for magnesium and of calcium for sodium, the dissociation of calcium carbonate in water, the effect of sulfate and undissociated calcium sulfate and magnesium sulfate on the exchange processes. Flow is assumed to be one-dimensional in the unsaturated zone and two-dimensional in the saturated zone. The model requires detailed input data, including:  (1) parameters that characterize the properties of both unsaturated and saturated zones; (2) boundary and initial conditions for flow; (3) concentrations of calcium, magnesium, sodium, chloride, bicarbonate, and sulfate in the spoil; (4) concentrations of the same species in water flowing into the spoil; (5) the cation exchange capacity; (6) quantity of gypsum; and (7) percentage of calcium carbonate present. Output consists of infiltration rate, volumes of subsurface drainage, and concentrations of calcium, sodium, magnesium, bicarbonate, chloride, sulfate, and dissolved solids in the drainage waters.
The model was tested using data from a site in southwestern Colorado. The model successfully predicted the infiltration rate andable
